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Abstract. 

The occurrence of a bounce in the FRW cosmology requires modifications of General 
Relativity. An example of such a modification is the recently proposed Hofava-Lifshitz theory 
of gravity, which includes a "dark radiation" term with a negative coefficient in the analog 
of the Friedmann equation. A modification of the HL gravity, relaxing the "detailed balance" 
condition, brings additional terms to the equations of motion, corresponding to stiff matter. This 
paper presents comparison of the phase structure of the original and modified Hof ava cosmology. 
Special attention is paid to the analysis of a wide range of bouncing solution, appearing in both 
versions of the Hof ava theory. 



1. Introduction 

There have been many attempts to modify Einstein's theory of gravity to avoid an initial 
singularity. Some were made at the classical level, some involve quantum effects. Examples 
include the ekpyrotic/cyclic model ([IJ [21 El HI El IE]) and loop quantum cosmology (0 HI IS])) 
which replace the Big Bang with a Big Bounce. Attempts to address these issues at the classical 
level include braneworld scenarios ( [10l HI] ) , where the universe goes from an era of accelerated 
collapse to an expanding era without any divergences nor singular behavior. There are also 
higher order gravitational theories and theories with scalar fields (see |12j for a review of bouncing 
cosmologies). However it is fair to say that the issue of the initial singularity still remains one of 
the key questions of the early Universe cosmology. 

Recently much effort has been devoted to studies of a proposal for a UV complete theory of 
gravity due to Hof ava [U2 [HI [15] and modifications of the theory [TU [17\ [TBI US] (f° r a recent 
review see [2D])- Because in the UV the theory possesses a fixed point with an anisotropic, Lifshitz 
scaling between time and space, this theory is referred to as the Hofava-Lifshitz gravity. From 
the time at which the Hofava theory was presented, there is also quite much discussion of possible 
problems and instabilities of Hofava-Lifshitz gravity [2U [221 [231 E3J- Numerous sophisticated 
versions contain new terms added to the original Lagrangian with attempt to make the proposal 
more general [TjJ] and to solve the so called strong coupling problem [2JJ [521 HS1 HZ1 12H] • Even so 
it is still tempting to investigate issues opened by this theory and its modifications. 

Soon after this theory was proposed many specific solutions of this theory have been found, 
including cosmological ones ([211 E01 Ell EH HS1 E3 EH ESI EH EI])- It was also realized that the 
analog of the Friedmann equation in the HL gravity contains a term which scales in the same 
way as dark radiation in braneworld scenarios 29, 30, 3 lj and gives a negative contribution to 
the energy density. Thus, at least in principle it is possible to obtain non-singular cosmological 
evolution within Hofava theory, as it was pointed out in [221 ED E21 SO] • Such possibility may 
have dramatic consequences for potential histories of the Universe - other than avoiding the 
initial singularity. New imaginable scenarios of cosmological evolution include contraction from 
the infinite size, bounce and then expansion to infinite size again, or eternal cycles of contraction, 
bounce and expansion. 
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Additional possibilities are brought by some interesting modifications modifications of HL 
gravity, either by softly breaking a detailed balance condition [T31 UHl [T7] or relaxing it completely 
[TBI [IS]. I n the latter works Sotiriou, Visser and Weinfurtner (SVW) in search for the more general 
renormalizable gravitational theory took a gravitational action containing terms non only up to 
quadratic in curvature, like in original HL formulation, but also cubic ones, as suggested earlier 
in |29] 130] . Generalized Friedmann equation of this model include modified dark radiation term 
proportional to ~ 1/a 4 (a is a scale factor) from the original HL formulation, and also additional 
1/a 6 term. This new term, negligible at large scales, becomes significant at small ones and modifies 
bounce solutions. Specially, with the opposite sign (the value of a coupling constant is arbitrary) 
than the 1/a 4 term, it may compensate the dark radiation term at small distances and cancel the 
possibility of avoiding singularity, like in the HL gravity with the softly broken detailed balance 
condition and negative spacial curvature |39j . Thus one of the questions to be answered is how 
the additional terms in the generalized Friedmann equations of the SVW HL gravity influence the 
existence and stability of a cosmological bounce. 

In this work we are going to investigate, with the help of the phase portrait techniques, how 
relaxing the detailed balance condition affects the dynamics of the system, and then compare the 
results to those in the standard HL theory. We will we focus on non flat cosmologies, with space 
curvature k = ±1, allowing non-singular solutions. Unlike in our previous paper [42] . we are going 
to describe matter by a cosmological stress-energy tensor added to the gravitational field equations. 
Such analysis is more effective and avoids unneeded approximations and simplifications. In this 
hydrodynamical approach two quantities: density p and pressure p describe matter properties. 

Nonetheless, constant parameter w of the equation of state is of course an idealization, hard 
to avoid at this level of research. It would be better to have history of the SVW HL Universe 
constructed in a similar way as in the standard ACDM model, with phases and epochs of different 
matter/radiation contents. Yet unless a rich structure of the original and the generalized HL 
theory, with additional coupling constants whose range of values is not fully understood thus far, 
is investigated deeper, we shall use simpler tools. Thus in this paper, within limited physical 
understanding of the theory and its parameters, we would rather present lists of possibilities than 
likely physical solutions. With the progress of research in this field and better understanding of 
the nature of these parameters, it will be possible to assign more physical interpretation to a set 
of solutions/scenarios found. 

Related analysis of the generalized Horava-Lifshitz cosmology have recently appeared in 43 
and [44] . which we become aware of while this work was being typed. Those papers address a 
somewhat different set of issues, i.e. static solutions of the HL universe. Here we are interested in 
stable and unstable solutions leading to cosmological bounce, and focus on both cases of non flat 
universe (k = — 1 and k = 1) with the value of a HL constant A arbitrary. We agree on the regions 
of overlap. General discussion of the full phase space of the original HL cosmology is contained 
e.g. in [4%ll46j . 

The structure of a paper is following: in Section 2. we briefly sketch the Horava-Lifshitz 
gravity and cosmology. In Section 3. the possibility of bounce in theory with detailed balance 
condition is discussed. In Section 4 we discuss phase portraits of the HL universe with condition 
of detailed balance relaxed. 

2. Horava-Lifshitz cosmology 

The metric of Horava-Lifshitz theory, due to anisotropy in UV, is written in the (3 + l)-dimensional 
ADM formalism: 

ds 2 = -N 2 dt 2 + g ij (dx i - N i dt)(dx j - N j dt), (1) 
where A, iVj and gij are dynamical variables. 
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2.1. Detailed balance 

The action of Hofava-Lifshitz theory is [2] 

I = J dtd 3 x(£ + £i), (2) 

-4A) 2 n 2 r 
32(1 -3A) 2u^ lJ 

where Kij — -h \_\fjij ~ V(i-^j)] i s extrinsic curvature of a space-like hypersurface with a fixed 
time, a dot denotes a derivative with respect to the time t and covariant derivatives are defined 
with respect to the spatial metric gy. Moreover 

pUJ 2 

Zij = Cij ^Rij- (3) 
K , A, p, w and are constant parameters and the Cotton tensor, C^, is defined by 

C ij = fl ki Vk f R j i _ ^RSj^j = e^VkR^ - ^e ikj d k R. (4) 

In ([2]), £ is the kinetic part of the action, while C\ gives the potential of the theory in the so-called 
"detailed-balance" form. 

Matter may be added by inserting a cosmological stress-energy tensor in gravitational field 
equation. Within such framework we approximate the stress-energy tensor by two quantities: 
density p and pressure p, then simply add them to the vacuum equations by demanding the 
correct limit as one approaches General Relativity - the low energy limit of the HL theory. Relation 
between p and p is given by the equation p = wp, with w being the equation of state parameter. 

Comparing the action of the Hofava-Lifshitz theory in the IR limit to the Einstein-Hilbert 
action of General Relativity, one can see that the speed of light c, Newton's constant G and the 
effective cosmological constant A are 

K 2 fi I A w n k 2 c 3n 4 p 2 A 2 w 3c 2 

C= ^VT^3A' G= 32V A = -3A^1^^^' (5) 
respectively. To have real value of speed of light c emerging, the HL cosmological constant has 
to be negative for A > 1/3 and positive for A < 1/3. It is possible to obtain a positive cosmological 
constant A^ in the IR limit A = 1 if one performs in ^ an analytic continuation of constant 
parameters p ip and uj 2 —iui 2 . 

The equations for Hofava-Lifshitz cosmology are obtained by imposing conditions of 
homogeneity and isotropy of the metric. The associated ansatz is N — N(t), 2Vj = 0, = a 2 (i)7^- 
where a(t) is a scale factor and 7^ is a maximally symmetric constant curvature metric, with a 
curvature k = {— 1,0,1}. On this background 
H 2k 

where H = a/a is the Hubble parameter. 
The gravitational action ([2| becomes: 



K-ij — AT 9ij i Rij — o^ij ' — 0, (6) 



f , ,3 „ 3 f 3 (! ~ 3A ) r2 ZK 2 p 2 k W ( k A w \ 

- "V #\ (7) 
8(l-3A)a 4 J K > 

The equations of motion are obtained by varying the action Q with respect to N, a and <p, setting 
N = 1 at the end of the calculations and adding terms with density p and pressure p, leading to 

2 = K 2 P K 4 p 2 ( kA W _ Ajjr _ _fc2_\ 

6(3A- 1) 8(3A- l) 2 V " 2 o„4 J ' 



kA\y 


A 2 h 2 

ft 


a 2 


~Y ~ 2a7 


( kA w 


k 2 \ 


• ( a 2 


+ 4^)< 



■ = n 2 (p+p) n 4 p 2 ^ .„ 

4(3A- 1) + 8(3A- l) 2 [ -2 1- „„ 4 J ' 
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and the continuity equation: 

p + 3H(p + p) = 0, 



(10) 



The upper sign denotes the Aw < case, the lower one the analytic continuation p ip with a 
positive Aw. 

The significant new terms in the above equations of motion are the (l/a 4 )-terms on the right- 
hand sides of (|8| and They are reminiscent of the dark radiation term in the braneworld 
cosmology |47] and are present only if the spatial curvature of the metric is non- vanishing. 

Values of constant parameters n 2 and p may be expressed in terms of cosmological constants 
according to We will also work in units such that 8ttG = 1 and c — 1. Then 



and 



32vrGc : 



±- 



A= -A w , 



1 - 3A 2A 

Substituting the above expressions and the equation of state p 



wp to leads to 



H 2 = 



H = 



3A- 1 
2 

3A- 1 



P 



± 



A 

3 ~ V 3 " 
(1 + w) 



3 


k 2 Y 




+ 4A 


oV_ 




(-- 


3 k 2 


\a 2 


2A~o7 



(11) 

(12) 

(13) 
(14) 



2.2. Beyond detailed balance 

The gravitational action written in the "detailed balance" form ^ ( 14J ) contains terms up to 
quadratic in the curvature. However the most general renormalizable theory contains also cubic 
terms, as it was pointed out in [29l [30]. Thus Sotiriou, Visser and Weinfurtner ( [THl [19]) built 
a theory with projectability condition N — N(t), as in original Hofava theory, but without the 
detailed balance condition. This led to Friedmann equations with an additional term ~ 1 /a 6 and 
uncoupled coefficients: 



H 2 



H 



(3A - 1) 



(3A - 1) 



■<Ti + tT2" 



P 

2 



P 
2 



(J2- 



cr 3 - 



2(7 3 - 



cr 4 - 



3(74- 



(15) 
(16) 



Values of constants 03 , 04 are arbitrary. In order to coincide with the Friedmann equations in the 
IR limit A = 1 and for large a (terms proportional to 1 /a 4 and to 1/a 6 are then neglible) one has 
to set cr 1 = A/3 and 01 = —1. Thus the above equations take the following forms: 



H 2 



H 



(3A - 1) 



(3A - 1) 



A 
3" 



k 



k 2 



0-3- 



(74- 



2(7 3 



3(74" 



(17) 
(18) 



where we have used the equation of state p — wp. Density parameter follows the standard evolution 
equation ( 10 ). New terms, proportional to 1/a 6 , appearing in the analogs of Friedmann equations, 
mimic stiff matter (e.g. such that p = p and P s yff ~ l/° 6 )- These terms are negligibly small at 
large scales, but may play a significant role at small values of a scale parameter. 



3. Bounce stability in the original HL theory 

In order to investigate the appearance of a bounce in the original HL gravity, we are going first to 
simplify the equations of motion (8|9| and then to reduce them with respect to Q. In this way 
we will obtain the two dimensional dynamical system describing the evolution of a and H. 
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Solving Eq. ( 13 ) for p gives 
3(3A 



P 



A-3- 



4A~aJ 



Inserting the above formula to (fl4| leads to 
±1 



H 



3A - 1 



x k 3(3w- 1) k 21 
(1 + id) A - (3w + 1) ^ + 



4A 



3 (1 
2 V 



w)H 2 . 



(19) 



(20) 



Equation (20) and the definition of the Hubble parameter: 
a = aH, 



(21) 

provide the two dimensional dynamical system for variables a and H . 

To find the finite critical points we set all right-hand-sides of equations ( 20pT ) to zero. This 
gives two points: 

' >U ' H = 0, (22) 



Pi : a = 



Pi 



2A' 
(3w- l)fc 



,H = 0. 



(23) 



(l+u>)2A' 

These points arc finite, unless w — — 1. In the latter case point P2 is moved to infinity. Point Pi 
exists for k/A > 0. Point P2 exists for w > 1/3 and k/A > or w < 1/3 and k/A < 0. Thus those 
two points exist both at the same time for w > 1/3. 

Stability properties of the critical points are determined by the eigenvalues of the Jacobian 
of the system ( 20||2"T |. More precisely, one has to linearize transformed equations ( 20||2~lj at each 
point. Inserting x = xq + Sx, where x — (a,H), and keeping terms up to 1st order in Sx leads 
to an evolution equation of the form Sx — ASx. Eigenvalues of A describe stability properties 
at the given point. Critical points at which all the eigenvalues have real parts different from 
zero are called hyperbolic. Among them one can distinguish sources (unstable) for positive real 
parts, saddle for real parts of different sign and sinks (stable) for negative real parts. If at least one 
eigenvalue has a zero real part (non-hyperbolic critical point) it is not possible to obtain conclusive 
information about the stability from just linearization and needs to resort to other tools like e.g. 
numerical simulation |48] . 

Eigenvalues at Pi are following: 





/ T2A \ 


/3(1-3A)' 2 V 


3(1 - 3A) J 



For all admitted values of A and A, expression =pA/(l — 3A) is negative, thus Pi is a center (both 
eigenvalues of A are purely imaginary at this point). 

Eigenvalues at the second finite critical point P2 read as: 



T2A(1 + w) 



,2 



T2A(1 + w) 



\ ]] (l-3A)(l-3to)' y (l-3A)(l-3tu) J' 

Depending on the value of parameter 10 the point P 2 may be a center or a saddle (two real numbers 
with opposite signs). Precisely, P2 is a linear center (non- hyperbolic center with one eigenvector) 
for id = — 1 (k/A < 0), a center for -1 < w < 1/3 (k/A < 0) and a saddle for w > 1/3 (k/A > 0). 
Properties of the critical points Pi and P2 in dependence on the values of A, k, w are gathered 



in the Table 1. A < corresponds to solutions of (20)-(|21| with the upper sign, the case of A > 



to the lower sign in ( 20 1-( 21 ) 



To find critical points that occur at infinite values of the parameters we rescale the infinite 
space (a,H) into a finite Poincare sphere (as in [SD]) m such a way that the new coordinates 
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w 


Pi 


stability 


Pi 


stability 


> 


>i_^ 


+ 


center 


+ 


saddle 


-1 < W < \ 


+ 


center 






-1 


+ 


center 






< 












-1 < w < k, 






+ 


center 


-1 






moves to oo 


linear center 



Table 1. Properties of finite critical points in the HL theory. The plus sign stands for "exists" 
and the minus sign stands for "does not exists" . 



(a, H) are written in polar coordinates r,<p: a — r cos <f> and H = r sin <j) and: 



r 



1-r 



COS0, (24) 



H=—^—sin<j), (25) 
1 — r 

We also rescale the time parameter t by denning the new time parameter T such that: dT = 
dt/(l — r). In these coordinates our phase space is contained within a sphere of radius one - 
infinity corresponds to r = 1. More precisely, semi-sphere, as a scale factor a may take only 
nonnegative values. 

This leads to the dynamical equations in terms of r, <j) and their derivatives with respect to 
new time T. Taking limit r — 1 we obtain: 

r'(T) = 0, (26) 
5 -4- ^7/) 

4>'{T) = - cos0sin 2 0. (27) 

Putting r.h.s. of the above equations to zero, we find 4 solutions: 

^3 = (1,0) 

P 4 = (1,tt/2) 

P> = (1,7T) 

P 6 = (1,3tt/2) 

in polar coordinates (r,4>). Point P5 is nonphysical (a negative a) and shall be removed from 
further discussions. Eigenvalues of the Jacobian matrix at the above points are following: 

(0, 0) at P 3 : 
5 + 3w 1 + ■„ . 
— 2~ ,3— 2~~ 
5 + 3w „ 1 



at P 6 

The point P3 is non-hyperbolic and we determine its properties by numerical simulations for each 
set of parameters. Unless w = — 1 points P4 and P 6 are respectively a repelling and an attracting 
node. For w = — 1 the finite fixed point P2 is moved to (00, 0) becoming P3, which is then a 
linear center. For this value of w points P4 and Pq are non-hyperbolic. It follows from numerical 
simulations that they are saddles then and ends of a separatrice. 

Fig. 1 shows the phase portrait of HL universe containing matter with equation of state 
parameter w > 1/3, k/A > 0, Figure 2. shows phase portrait for — 1 < w < 1/3, k/A < and 
Fig. 3 for w = — 1, k/A < 0. One has to note that these figures contain the deformed phase space, 
scaled to fit on the finite Poincare sphere. One may have the impression that they describe regions 
in which e.g. the scale factor a increases although the Hubble parameter H is negative. However 
it is the parameter a that is increasing on the diagram, not the scale factor a. 

Bounce scenarios are thus possible when critical points exists. If these point are centers then 
there are closed orbits around them and the Universe goes through eternal oscillations: expansion, 
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Figure 1. Projected phase space of HL universe with kA > and w > 1/3. 
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a 

Figure 2. Projected phase space of HL universe with k/A < and —1 < w < 1/3. 
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collapse to a finite size, expansion etc. Point P± may be a center (for certain values of parameters), 
but then p — 0, which is physically not interesting. More interesting case is when P 2 is a center, 
then there are closed orbits with a non-zero density p. The third bounce scenario is around the 
linear center P2 (moved to 00 and coinciding with P3). In this case there are onefold closed 
trajectories: universe starts from a static one (H = 0) and infinite (a = 00), goes through a period 
of collapsing to a finite size, then after a bounce starts expansion and finishes as a static infinite 
universe. 



4. Bounce stability in the SVW generalization 



Relaxing detailed balance condition leads to the generalized Friedmann equations ( 17p8 ) with 
additional term ~ 1/a 6 and uncoupled coefficients. 
We may solve eq. |l7| ) for p obtaining: 

„(3A-1) 2 k a 3 k 2 3a 4 k 

p = 3 H -A-3^-3 — 5 -g-. 



(28) 



Substituting this expression on p into (18) and using the equation of state p = wp leads to 
2 M(l + w) fe(l + 3«;) 
~~2 2a 2 



i? = 



3(1 + w) h2 
2 



(29) 



3A - 1 

cr 3 (-l + 3w)fc 2 3a 4 (l + w)fc x 
+ 2a 4 _+ 2oT 

The above equation, together with the definition of the Hubble parameter provides the two 
dimensional dynamical system for variables a and H . 

Finite critical points are solutions of the equations (21 1 and (29 1 with r.h.s. set to zero. Hence 
these points fulfill H = and: 

A(l + w)a 6 - k{\ + 3w)a i + cr 3 (-l + 3w)k 2 a 2 + 3ct 4 (-1 + w)k = (30) 
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4-1- Cosmological constant w = — 1 



For w = — 1 equation (30) reduces to a biquadratic equation: 
ka — 2osk 2 a 2 — 3a^k = 0. 



(31) 



Solutions of (31) are following: 
Pi : a 2 = fc(T3 — 
P, : a 2 = fcer 



: + 3a 4 , (32) 

■2 '■ il' = l'' n ':; - \l &3 +3(74. (33) 

Point Pi exists when {(ko-3 > 0, 04 < 0); (|ct4| < er|/3)}. Point P2 exists for: 

{(kcr 3 > 0,cr 4 > 0); (fccr 3 > 0,cr 4 < 0, 1 0-4 1 < of/3); (fca 3 < 0, cr 4 > 0)} . (34) 
Stability properties of the critical points found are given by eigenvalues of the Jacobian matrix A 



of the system (|2l|),(|29|). Eigenvalues of A at Pi are following: 

-kCi _ / -jfcCi \ 
3A^T ' 




(35) 




where Ci denotes expression in fc,cr3,(74, positive when the point Pi exists. Similarly, eigenvalues 
of A at the point Pi read as: 



(36) 



where C2 denotes expression in k, 0-3,0-4 being positive when the point P2 exist. 

Thus for fc/(3A— 1) > the point Pi is a center and P2 an unstable saddle, for fc/(3A— 1) < 
Pi is a saddle and P2 a center - provided that the values of k, 03, 04 allow their physical existence 
(a 2 > 0). 

Density p at those points is equal to 



P = 



2o-f + 9er 3 cr 3 + 9Laj ± 2fc(j|^cr| + 3f7 4 ± 6fe<T4-v/o| + 3f7 4 



9<rl 



(37) 



where '+' corresponds to the point Pi and '-' to P 2 . Thus at Pi density p > if this point exists 
and (k < 0, cr 3 < 0, > A > l/er 3 ). At the point P2 density is positive if (k < 0, a 3 < 0, A < 0). 

^.12. Radiation w = 1/3 



When ui = 1/3 equation (30) reduces to the following one: 
2A 



fcx — koA 



0, 



(38) 



where x - 
form: A 
as: 



The discriminant of the cubic polynomial a 3 y 3 4- Q2y 2 + aiy + ao is of the following 



18a aia2a 3 — 40^0 + a 2 ; ! ~" 4a 3 af — 27agag. Discriminant of the equation (38) reads 



A = -4fc 2 t7 4 (fc 2 + 3A 2 o- 4 ) 



(39) 



For A > the cubic equation (38) has three real solutions. Condition A > is fulfilled for a 
nonflat universe (k 7^ 0) when 04 < and \<ta\ < 1/(3A 2 ). Otherwise (38) - the equation with real 



coefficients - has one real solution and two nonreal complex conjugate roots (A < 0) or multiple 
real roots (A = 0). 
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Physical points exist if real solutions x = a > 0. Equation ( 38 ) cannot have three real positive 
roots, as it is implied by Viete's formulas - precisely, by the second formula of the following ones: 

3fc 

— = Xi + x 2 + x 3 , (40) 
= X1X2 + X2X3 + X3X1, (41) 

3fc(T4 

-^=x 1 x 2 x 3 . (42) 

If there are three real solutions (A > 0), one or two of them may be also positive. There is 
one positive solution if k/A < and two positive solution if k/A > 0. 

Multiple real solutions exist when A = hence when 174 = or 0-4 = — fc 2 /(3A 2 ). The 
former case corresponds to HL theory with the detailed balance condition, the latter case implies 
solutions: 

Qx : a 2 = -k/(2A), (43) 
Q 2 : a 2 = k/A (double root), (44) 

and H = 0. Depending on the sign of k/A one of the two solutions has physical meaning. 

One real solution (A < 0), i.e. when cr 4 > or o~4 < — 1/(3A 2 ), has positive value if ka 4 /A > 0. 
Eigenvalues of the Jacobian matrix at the critical points a x are following: 



2 k(3a 4 + ai) 2 /fc(3cr 4 + <4) 



3A - 1 ' al V 3A - 1 



(45) 



Thus they may be saddles or centers, depending on the sign of k/(3\ — 1) and 3o~4 + a*. 

We can describe more precisely the case when 0-4 = — fc 2 /(3A 2 ) (0-4 = case is described 
within original HL cosmology) and critical points are a 2 = — fc/(2A) (Qi) or a 2 = k/A (Q2)- Then 
the eigenvalues of Jacobian matrix read as: 



(WaA^i'-WaA^ij at (46) 

and 

(0,0) at Q 2 . (47) 

Therefore the point Q\ may be a saddle or a center, depending on the sign of A/(3A — 1). The 
points Qi is non- hyperbolic, numerical simulations (Fig. 5) show that it is a cusp. Comparing the 
neighborhood of Q2 on both this figures, one can see the difference between the deformed phase 
space and the non-deformed one. At the former one there are regions in which the parameter a 
increases although the parameter H is negative, whereas this behavior is absent at the latter one. 
This is due to the fact that parameters a and H are geometric objects without the same physical 
meaning as the scale factor a and the Hubble parameter H. 
Density p is following: 

p= -3A(5 + 4Acr 3 ) at Qi, (48) 

p= -3A(-1 + Act 3 ) at Q 2 . (49) 

Density at Qi is positive if 

{(A < 0, 0-3 < 0); (A > 0, 5/(4A) > a 3 )}. 

At Q2, p is positive for 

{(A<0,Aa 3 >l);(A>0, ( T3<l/A), 
plus conditions for existence: Q\ has physical meaning if k/A < and Q 2 if k/A > 0. 
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4-3. General case 

In general critical points of the system (21) and (291 are of the following form: (a x ,0), where a 2 



is a root of the equation: 

A(l + w)x a - k(l + 3w)x 2 + (7 3 (-1 + 3w)k 2 x + 3(7 4 (-l + w)k = 0. (50) 

Depending on the sign of the discriminant A = 18(10010203 — Aa^aQ+a^al — 4030^ — 27a 2 ag, where 
ao = A(l + w), ai = — fc(l + 3iw), 02 = os(— 1 + 3w)k 2 , 02 = 2>o-^{— 1 + w)fc, the above equation 
has one, two or three real solutions. Namely, for A > 0, there are three real roots, for A < there 
is one real root and two complex conjugates, for A = those conjugates become a real double 
root. 

The eigenvalues of Jacobian matrix at those critical points read as: 



2Difc / 2D x k 



o6(3A-l)'Vo6(3A 



(51) 



where 

D 1 = (l + 3w)a 4 + 2fccr 3 (l - 3w)a 2 - 9a 4 {w - 1). (52) 

For w < —1/3 and —(74 > cr|(l — 3w) 2 /((3w + l)(w — 1) expression D\ is always negative and it 
is always positive for w > —1/3 and 0-4 > cr|(l — 3w) 2 /((3w + 1)(1 — w). Thus depending on the 
sign of fc/(3A — 1) the critical points, if exist, are either always stable or always unstable. Nature 
of their stability depends on the values of a x , A, a 3 and cr 4 . 

Stability properties of critical points at infinity is the same as for the detailed balance case. 



After the Poincare transformation |24|) and ( 25 ) the whole phase space is contained within a 
semi-circle (a > 0) of radius one. 

Points at r = 1 and cf> — tt/2, 3tt/2 are repelling and attracting node, respectively. Point 
r = 1, (f> = is non-hyperbolic, and its stability properties can be obtained e.g. from numerical 
simulations. 

Figure 6. shows the phase space of system with three finite critical points. Points S± and S3 
are centers, point S2 is a saddle. 

In Table 2. we have gathered properties of the finite critical points in the SVW generalization 
of the Hofava cosmology. 



5. Conclusions 

In this work we have performed a detailed analysis of a phase structure of the HL cosmology with 
and without detailed balance condition. Both this models contain a dark radiation term 1/a 4 
in the analogs of the Friedmann equations. Thus it is possible for a nonflat universe (k ^ 0) 
that the Hubble parameter H = at some moment of time, which is a necessary condition for 
the realization of the bounce. Comparing phase trajectories obtained in those models we have 
attempted to answer the question how the generalization of Hofava gravity (breaking the detailed 
balance condition) impacts the occurrence and behavior of bouncing solutions. Additional term 
1/a 6 that appears in the Friedmann equations of SVW model, is of either sign, and thus it may 
possibly compensate the 1/a 4 term (generic for HL gravity) leading to the singular solution. 

Indeed, it occurred that the biggest difference between the Hofava theory and its 
generalization arrives for the small values of a scale parameter a and a Hubble parameter H. 
This is not surprising, as the SVW gravity term 1/a 6 plays role only for the small values of a and 
becomes insignificant for the bigger ones. 

In the original Hofava formulation there may be two finite critical points, one of them a 
center and one a saddle. They are of the type (a^O) in the (a, if) space, thus strictly connected 
to bounce solutions. These pairs of points exist both only for matter with w > 1/3. Around a 
center there are closed orbits corresponding to the oscillating universe, i.e. going through eternal 
cycles of contraction, bounce and expansion. These orbits resemble bounce solution described by 
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0.0 0.2 0.4 0.6 O.S 1.0 



Figure 6. Projected phase space of the generalized HL universe with 3 critical points existing 



Brandenberger [55] or quasi-stationary solutions presented in our previous work [32] based on a 
field approach. Such solutions are physically interesting (density p > 0) for w < 1/3 and k/A < 
- so either for a closed universe with a positive cosmological constant, or an open universe with 
k = 1 and a negative cosmological constant A. The second class of oscillating solutions, with 
vanishing density p = 0, appears when k/A > 0. Additionally, there is a third bounce scenario 
around a linear center P 2 (moved to oo and coinciding with P 3 ). Here there are onefold closed 
trajectories: universe starts from a static one (H = 0) and infinite (a = oo), goes through a period 
of contraction to a finite size, then after a bounce starts expanding and again ends as a static 
infinite universe. Moreover, for some values of parameters, i.e. k/A < and w > 1/3 there are no 
finite critical points, thus no bouncing solutions. 

In the SVW HL cosmology, with additional term appearing in the analogs of Friedmann 
equations, there may exist 0,1,2 or 3 finite critical points. They are also of the type (a x ,0) in 
(a,H) space. Here a x is the solution of the bicubic equation. In general there exists at least 
one real solution of the cubic equation with real coefficients, but physical points correspond only 
to positive values of these roots. Critical points might be stable centers - surrounded by closed 
orbits, describing oscillating universes, or unstable saddles. There also exist solutions with orbits 
around a linear center at (oo,0), where similarly as in the original HL theory, a universe starts 
from a static infinite the, collapses to a finite size, undergoes a bounce and then expands to a 
static infinite state. Thus there is one cycle only, without further oscillations. There are also sets 
of parameters, much wider than in the original HL theory, that do not allow the existence of finite 
critical points, leading only to singular solutions. 

The most significant feature of oscillating (and bouncing) solutions in the SVW formulation 
is the existence of two centers, with a saddle between them (three finite critical points) for some 
values of parameters. We present such a solution at the Fig. 6. In a more realistic situation, 
that includes dynamical change of state parameter, it would be possible to go from one oscillating 
bouncing solution to another. Present framework does not allow such evolution as it describes 
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Point 


w 


Existence 


Stability 


p positive 


Pi 


-1 


(ka 3 > 0,0-4 < 0) 


center for 
fc/(3A - 1) > 

saddle for 
fc/(3A - 1) < 


k < 0, 0-3 < 
> A > 1/0-3 


P2 


-1 


(fccr 3 > 0, 
> cr 4 > -0-3/3) 
{ka 3 < 0, er 4 > 0) 


center for 
fc/(3A - 1) < 

saddle for 
fc/(3A - 1) > 


k < 0, cr 3 < 
A < 


Qi 


1/3 


((74 = -fc 2 /3A 2 , 
fc/A < 0) 


center for 

A/(3A- 1) < 

saddle for 
a /(*> \ 1 \ ^ n 


(A < 0, (73 < 0) 
(A > 0,0-3 <-5/(4A)) 


Q2 


1/3 


((7 4 = -fe73A 2 , 
fc/A > 0) 


cusp 


(A < 0,Acr 3 > 1) 
(A > 0,(7 3 < 1/A) 


Si 
S 2 


(-1,-1/3) 


depending on 
the sign of 

/\ fmn at i"Vir» rn f~\T c • 
Z_l cLlJAl Ul lUUbo. 

0,1,2 or 3 points 


" 4 ^ (3w + l)(w-l) 

center for 

saddle for 
fe/(3A-l) < 




Si 

S2 
S3 


(-1/3,1/3) 
U(l/3,1]) 


depending on 
the sign of 
A and of the roots: 
0,1,2 or 3 points 


a4 > (3w+l)(l-w) 

center for 
fc/(3A- 1) < 

saddle for 
fc/(3A - 1) > 





Table 2. Properties of the finite critical points in the SVW HL theory. 



matter as hydrodynamical fluid with a constant w. We expect that the field approach, with a 
more complete dynamics, may be suitable for further investigation of this interesting scenario. 

The phase structure at infinity is the same at both formulations. Except bouncing solutions 
around finite critical points, there are also solutions leading to Big Bang, Big Crunch or eternal 
expansion. It is worth to stress that in both models, the original HL gravity and the SVW 
generalization, there are classes of parameters that do not allow a non singular evolution. Physical 
interpretation of some of these parameters (coupling constants 03 and (74 in SVW model) still 
remains an open question. 
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